Helioseismic Frechet Traveltime Kernels in Spherical Coordinates 

R. B. Schlottmann and A. G. Kosovichev 
W. W. Hansen Experimental Physics Laboratory, Stanford University 

Stanford, CA 94305 

3 ' briansch@f usemail . com 

^ ■ Received ; accepted 



43 

9 L,: 
6 

(N 
> 

O 



To be submitted to The Astrophysical Journal 



-2 - 
ABSTRACT 

Seismic traveltime measurements are a crucial tool in the investigation of 
the solar interior, particularly in the examination of fine-scale structure. Tra- 
ditional analysis of traveltimes relies on a geometrical ray picture of acoustic 
wave propagation, which assumes high frequencies. However, it is well-known 
that traveltimes obtained from finite-frequency waves are sensitive to variations 
of medium parameters in a wide Fresnel zone around the ray path. To address 
this problem, Frechet traveltime sensitivity kernels have previously been devel- 
oped. These kernels use a more realistic approximation of the wave propagation 
to obtain a linear relationship between traveltimes and variations in medium 
parameters. Frechet kernels take into account the actual frequency content of 
the measured waves and, thus, reproduce the Fresnel zone. Kernel theory has 
been well-developed in previous work on plane-parallel models of the Sun for 
use in local helioseismology. Our primary purpose is to apply kernel theory to 
much larger scales and in a spherical geometry. We also present kernel theory 
in a different way, using basic functional analytic methods, in the hope that this 
approach provides an even clearer understanding of the theory, as well as a set 
of tools for calculating kernels. Our results are very general and can be used 
to develop kernels for sensitivity to sound speed, density, magnetic fields, fluid 
flows, and any other medium parameter which can affect wave propagation. 

1. INTRODUCTION 

In both the Sun and the Earth, the passage of waves, acoustic or elastic, respectively, 
through the body yields an indispensable probe of the inner structure and physics that 
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would otherwise be unavailable. Both helioseismology and geoseismology have relied on the 
inversion of normal mode data and time-distance measurements to obtain information on 
their subjects' interiors. In both fields, the time-distance measurements are the data that 
are more sensitive to small-scale, three-dimensional variations in structural parameters such 
as wave speed and density. 

Traditionally, the measurements of traveltimes have been interpreted using a 
geometrical ray approximation to the wave propagation. However, it has long been 
appreciated that this approximation is generally inadequate, geometrical rays being 
accurate representations of waves only at very small wavelengths relative to variations 
in medium parameters. In reality, the arrival time of a wave packet associated with a 
particular ray path is sensitive to structure well off the ray, the width of that "Fresnel zone" 
being dependent on the dominant frequencies in the wave. 

In an effort to obtain better information on the Earth's interior , geoseismologists 



Dahlen et al. 



2000 



Hung et al. 



20001) 



developed traveltime sensitivity kernels (See, e.g., 
which account for the finite frequency nature of the waves. By more accurately representing 
the physics of wave propagation, these kernels provide a better "map" between structure 
and traveltimes, which, in turn, yields a more accurate inversion result. 



Traveltime kernels are not new to helioseismology. 



Gizon fc Birch 



fl2002h developed 



2D kernels f or f-m ode traveltime sensitivity to acoustic damping and excitation, while 



Birch et al. 



(120041 ) . using a plane-parallel model of the Sun's near-surface, developed 
3D traveltime kernels for perturbations in sound speed and density. In both papers the 
theoretical foundations of sensitivity kernels, as well as practical considerations in their 
calculation, were laid out clearly and accurately. We will rely most heavily on those two 
sources to guide us in this work. 



The primary purpose of this paper is to further generalize the work of 



Gizon fc Birch 
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( 120021 ) and iBirch et al.l ( 120041 ) to a spherical geometry, whi ch would al 



ow th e development 



20091 ) . However, we 



of kernels for deep Sun structure, an area of active interest (jZhao et al. 
also wish to present the theory using the tools of basic functional analysis, with the hope 
that this way of looking at kernels provides an even clearer understanding of their nature. 

Our theoretical results are quite general and can be used to develop kernels for various 
structural parameters. Not only can formulas for sensitivity to sound speed and density be 
obtained, but one can also derive kernels for magnetic fields, fluid flows, and damping. In 
principle, the sensitivity to any parameter which can affect the propagation of waves can be 
derived. 



2. INTRODUCTION TO KERNEL THEORY 

In this section, we lay out a fairly general development of helioseismic traveltime 
sensitivity kernels based on the Born approximation. In geoseismology, where their 
theory and application were first developed, they are frequently called Frechet kernels or, 
more colorfully, "banana-doughnut" kernels, because of their characteristic appearance in 
cross-section. As will be made precise later, these kernels quantify the first-order (i.e., 
linear) dependence of seismic traveltimes on perturbations of any physical properties of a 
medium that can affect the propagation of seismic waves. By perturbation, we mean the 
difference between the value of a property for a particular medium and that of a reference 
medium. Usually, the reference medium is taken to be one with a high degree of symmetry, 
which greatly reduces the computational cost of calculating kernels. In the case of both the 
Sun and the Earth, this means that the reference model is most often spherically symmetric. 

The theoretical results contained in this section are not particularly novel. Kernels have 
been developed previously for the Sun in Cartesian coordinates and for the Earth (which 
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has somewhat different wave propagation effects from the Sun) in spherical coordinates. 
However, the development presented here utilizes a different theoretical apparatus which 
may be clearer and cleaner in some ways than the usual derivations. In particular, heavy use 
is made of some basic concepts in functional analysis. These concepts are not complicated, 
but since they may not be familiar to eve ry reader, we give below a non-rigorous review of 



what will be needed. (See 



Hatfield 



1992 



ch. 9, for a nice review of functional calculus.] 



2.1. Basic Functional Analysis 

First, we define a functional as a "function over function space" or, more precisely, a 
mapping from a space of functions to a finite-dimensional space of real or complex numbers. 
If / is a functional of a function g and itself a function of a variable x, we denote this as 
f[g](x). A relevant example would be the acoustic displacement wavefield, u, in a medium 
of sound speed c. Since the field is a function of position and time as well as a functional 
over the space of sound speed models, we would denote it here by u[c](x, t). 

The workhorse in our use of functional analysis will be the functional derivative, which 
is formally (i.e., non- rigorously) defined as 

SJM _ hm f[9(-)+eS(--a)]-f[g(-)] 

5g(a) e ->o e ' 1 1 

where the dot in g(-), for example, is a placeholder for the unspecified argument of function 
g. It may be useful to conceptualize the functional derivative as a gradient in an infinite 
dimensional space, with the value of the function at each point in physical space viewed as 
one coordinate of the function's "position" in function space. 

One consequence of our definition is that if 

/"OO 

f[g}(x) = / dx'g{x')8 (x - x')) = g(x), (2) 
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then our definition of the functional derivative implies 

Sf[g](x) 6g(x) 



89(a) 

which is a result we will use often. 



89(a) 



5 (x — a) 



(3) 



The only other concept from functional analysis we need is that of the functional 
Taylor series. If / is a functional over some function space, then / can be expanded about 
some specific function g : 



fig] = f[go 



where 



dx'Ag(x 



„ 8f[g] 



8g(x') 



9o 



-I /*0O /*00 

- / dx" / dx'Ag{x')Ag{x") 
^- J —00 J —00 



8 2 f[g] 



5g(x")5g(x') 



+ 



Ag(.) = g(.)-g (.). 



(4) 
(5) 



Note that all the above concepts relationships extend in an obvious way to functionals 
over functions of several variables as well as to functionals over vector functions. 



2.2. Development of the Kernels 

Here we define the inverse problem to be solved and present the mathematical 
derivation of the traveltime kernels that arise from it. We first cover the definition of 
some basic observables, the precise definition of what we mean by "traveltime", and the 
linearization of the relationship between traveltimes and actual medium properties. We 
then connect the results explicitly to wavefields in a reference solar model, providing general 
formulas for the kernels in terms of observational constraints (such as filters, line-of-sight 
effects, etc.), the reference wavefields, and the differential operators that govern them. 
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2.2.1. Some Definitions and the Linearized Problem 



Some care is required here in our definition of various quantities of interest since 
there are different observables which could be called the "wavefield" and different derived 
quantities which could be called the "data". We will refer to the actual particle motion 
due to the propagation of helioseismic waves in the Sun as the "wavefield" . The measured 
line-of-sight dopplergrams will be called the "raw data" or "raw traces", whereas the 
seismograms acquired by cross-correlating raw traces will be denoted "processed traces" or 
"processed data". 

In addition, we make a couple of brief comments about notation. To promote clarity, 
we indicate solar surface coordinates and integration variables by cr, cr', and so forth, 
whereas interior coordinates and integration variables are denoted by x and its variations. 
Also, unless otherwise stated, the Einstein summation convention is implied by repeated 
Latin indices. 

We denote the particle velocity wavefield by v(x, t) and the raw data by ip(cr, t). These 
two are related by 



where F represents the action of all spatial and temporal filters, which are assumed to be 
linear and time-translation invariant, that result from the instrument response and from 
whatever other filters that may be applied later. Note that the direction of F is along the 
line of sight from the point cr' to the instrument and that the integration is over the surface 
of the Sun. 

The processed traces derived from the raw traces are denoted D(cr a ,cr b ,t) and are 
obtained by cross-correlation of the raw traces associated with points er a and a b on the 
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solar surface: 



er a ,o- bj t)= dt'ip(er a ,t + t')ip((T b ,t') 



(7) 



J —oo 



For simplicity, we will not account for the finite time length of raw traces but instead will 
approximate the cross- correlation as an integral over all time. Note that we have defined 
D(cr a ,tTb,t) in such a way that for positive times, the resulting seismograms represent 
waves travelling from cr b to a a . 

We also define the quantities vo, ipo, and -Do a s the wavefield, raw data, and processed 
data, respectively, generated synthetically from a given reference solar model. Hereafter, 
any quantity with a subscript is taken to be associated with the reference model. 

We next define C(cr a , cr b ,t) as the windowed cross-correlation between the real and 
reference processed traces: 



C(a a , <r b , t) = / dt'W(tr a , tr b , t + t')D Q {a a , tr b , t + t')W(tr a , tr b , t')D(<T a , tr b , t'). (8) 



The function If is a windowing function applied to both the real and reference data, the 
purpose of which is to isolate pulses associated with specific propagation paths. 

Now let r(cr a , <j b ) be the time lag at which the cross-correlation above is maximum. 
We define this time as the traveltime. It is imperative to note that the traveltime must be 
measured this way before the kernels we derive here can be used. Other ways of measuring 
the traveltime are not necessarily better or worse, but they will require different kernels. We 
also point out that the way we have defined the cross-correlation ensures that if a pulse in 
the real data is early with respect to the reference data, then the traveltime will be positive. 
In other words, r is a measure of how much the real pulse is advanced with respect to the 
reference. 

The traveltime is a functional of the actual solar model. If we let q and qo denote 
vectors of functions describing the actual and reference solar models (e.g. 9i(x) and (/2( x ) 
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could be, respectively, the sound speed and density as a function of space inside the real 
Sun), then we can expand r in a functional Taylor series about qo, 



f St 
T[q](cr a ,(T h ) = T[q ]((T a ,cr b ) + ^ / d 3 x Ag Q (x) - — — (<r a , cr 6 ) 

„ J& oq a [x.) 



i V / d 3 x/ d 3 x'A ?Q (x)Ag^(x') — -, — . 
2! t^o J& <%(x')<9g Q (x 



qo 



+ 

qo 



(9) 

where Ag a (x) = q a (x) — g ,a( x )) an d integration is over the volume of the Sun. The first 
term in the expansion is identically zero since the time lag between the real and reference 
data would be zero if the real model were the same as the reference. In order to get a linear 
inverse problem to be solved later, we need a linear relationship between the traveltime 
and the perturbations Aq a to the solar model. This is accomplished by simply dropping 
the nonlinear terms in the above Taylor series. This is the only linearization required in 
this development of the kernels. Any other linearizations that occur later (e.g. the use 
of the first Born approximation) are a direct consequence of this one approximation and 
do not represent further approximations. However, this linearization does assume that 
the model perturbations are small in the sense that their main effect is to cause small 
changes in the arrival times of pulses in the real data. Larger perturbations invalidate this 
assumption because they can cause distortions in waveforms, not just time shifts, that make 
the cross-correlation traveltime measurements meaningless. Thus, hereafter we assume the 
equivalence 



r(cr a ,(T b ) = } d 3 x. Aq a (x) - — — (er a ,cr b ) 



(10) 



qo 



Traditionally, Frechet kernels are defined in such a way that they provide a linear 
relationship between traveltimes and fractional perturbations in model parameters. 
Following this tradition, we define K a , the kernel for the a-th model parameter, as 



8t 

K a (x; er a , a b ) = g ,a(x) t — rr^, cr b ] 

bq a {x) 
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qo 
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so that 



r(<T a , <r b ) = J2 I rf 3 x ^^(x; <r a , tr b ). (12) 



The remainder of this section will be devoted to relating the above expression, eq. (TTTj) . to 
calculable physical quantities and known observational constraints. 



2.2.2. General Kernel Formulas 



The expression in eq. fTTTT) is useless until we connect it explicitly to the underlying 
physics governing helioseismic wave propagation. We begin by providing an explicit 
mathematical definition of r. We have previously defined it as the value of the time lag that 
maximizes C, the windowed cross-correlation between the real and reference data. Letting 
C{t) = C(cr a , cr b , t), we expand C in a Taylor series about t = 0: 



t" 



(13) 



n=0 



Making the reasonable assumption that C is a smooth function, we have that C'(r) = 0. 
Thus, 



T n-l 



n=l 



n-1 ! 



C (n) (0) = 0. 



(14) 



Taking the first functional derivative of eq. (j!4p with respect to q a at q = qo, we get 







6C'(t) 



5C'(0) 
qo %*( x ) 



E 



-n-2 



Qo n=2 

Since r[qo] = 0, the above reduces to 

8C'(0) 



(n-2) 



■C (n) (0) 



5r 



r n-l ^(n)( ) 

5g a (x) (n - 1)! 5g a (x) 



(15) 



qo 







5g a (x) 



57 



+ Cff(0) 

q ^«( X ) 



qo 



which implies 



5g a (x) 



<5C'(0) 



qo 



cgf(o) 



(16) 



(17) 



qo 



where 



/oo 
dt'W(a a , a b , t + t')D ((T a , a b , t + t')W(a a , a b , t') A)K, tr b , t'). (18) 
-oo 

Eq. (|17p connects our definition of the kernel in eq. ffTTT) to something readily related to 
actual wave propagation. 

At this point, it will be more convenient to work in the frequency domain. To begin, 
we state our Fourier transform convention. If f(t) is a function of time, then 

1 



2tt 
1 



dte lut f(t) 



oo 
oo 



due~ luJt f(uj 



(19a) 
(19b) 



Now, if we define 



U n (cr a ,cr b , t') 



dt r 



W(tT a , CT b , t + t')D {cT a , (T b , t + t') W{(T a , (T b , f) 

L J t=o 

where n = 1,2, then from eq. (jHJ) 

/oo 
dt' ^i(ff o ,<T 6 ,0^K, 0-6,0 
-oo 

/oo 
du U*(cr a , a- bl u)D(cr a , <T b ,u), 
-oo 

and, similarly, 



(20) 



Co(0) 



du) U%(<T a , <T b , U})D (tT a , <T b , U)j, 



(21a) 
(21b) 

(22) 



where we have used the fact that and D are real-valued. With eq. ( I21b|) in hand, we 
find 



6C'{0) 



<5g a (x) 



qo 



duo U*(cr a , er b ,u] 



5D((T a ,(T bj Uj) 



5q a (x) 



(23) 



qo 



Thus, combining eqs. ffTTj) . fTTTl) . and f[2"3"j) . we have as an intermediate result 



if a (x; a a , a 



a, 



' cgf(o) 



do; i7*(<T , er b ,u/ 



5D(o- a ,cr b ,uj) 



(24) 



qo 



Since we will deal only with linearized wave propagation, we know the wavefield obeys 



£T(x, to) = S(x, oj) 



(25) 



where L is some linear differential operator that depends on the model parameters q, S is a 
source function, and T is a four-dimensional vector comprised of the particle displacement 
wavefield u and the pressure perturbation P' . If G(x, w;x ) is the Green's tensor in the 
actual solar model, i.e., if 

£G(x, cu; Xo) = S (x — x ) I, (26) 
where I is the identity matrix, then the solution to eq. (|25p can be written as 



T(x,w)= / d 3 x / G(x,w;x / )-S(x / ,w). 



(27) 



With this, eq. (Ej) can be rewritten in the frequency domain as 



4j{ct,u) = V27T / dV \ d 3 x.'F i (a-,(r',u})G ij {a',u}]x.')S j {x.',u) 

JdO JQ 



(28) 



where we take F to have zero pressure component. Using this in the definition of D, eq. ((?]), 

we get 



d cr\ I d xi / d cr 2 d x 2 F i (cr a ,cr 1 ,uj)Gi j (cr 1 ,u;x. 1 )S j (x.i,u) 

JO JdQ JQ 

x F£((T b ,(T 2 ,uj)Gl l ((T 2 ,uj;x 2 )S?(x 2 ,uj). (29) 



Regarding the source terms appearing above, we note that it is understood that the 
acoustic excitation in the Sun is due to turbule nt convection just below the photosphere 



and are often modeled stochastically. Following 



Birch et al. 



(120041 ) . we choose to represent 



the excitation by a source covariance matrix, M, defined as 

/oo 
dt'Sifct + t^Sjix'^t' 
-oo 



(30) 
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or, in the frequency domain, 



Mij(xL,x',u)) = Si(x,u)S*(x',uj). 



(31) 



(The value of using this representation will come later when we use a fairly simple reference 
model of the ensemble average of M. However, we will also consider the possibility that the 
true form of M could be taken as one of the model parameters in q.) Using this definition, 
eq. (1291) becomes 

D(cr a ,a b ,u) = 



d cr 1 cfxi / d er 2 d s x 2 Fi(o- aj cr ll uj)G ij (er l ,uj;]t 1 ) 

(90 J& J dQ JQ 

X F fc *(<T b , (T 2 ■ u)G* kl (tT 2 , w; x 2 )M,-j(xi, x 2 , u) 



(32) 



Within eq. (I32p . only the Green's functions and the source covariance matrix are ever 
dependent on the model parameters. Thus, for simplicity, when we take the functional 
derivative of D, we separate those out and define the quantity 

S 



d 3 xi / rf 3 x 2 Gi j (cri,a;;xi)G^(cr2,w;x2)M J 7(xi,X2,a;; 
<$g a (x) J Q J Q 



qo 

(33) 



In order to evaluate this, we will need the functional derivative of G. We start by taking 
the derivative of both sides of eq. ( 126]) , bearing in mind that the right hand side does not 
depend on the model parameters, 



5C 



<5G(x',u;;xo) 



G (x',u;;xo) + C — -—^ 

q ^aW 



0. 



(34) 



qo 



where Co and Go are the wave equation differential operator and Green's function, 
respectively, for the reference model. The solution to this equation is 



5G(x',u;;xo) 



Sq a (x) 



d 3 x"G (x',w;x") 



qo 



o 



6C(x",cu) 



Sq a (x) 



G (x",u;;xo) 
qo 



(35) 



-14- 



Putting all of our results together so far, we quote our general equations for the desired 
Frechet traveltime kernels: 



if Q (x; cr a , er b ) 



rvlt q °' a X n x dco d 2 cr 1 / d 2 a 2 Ui(cr a ,a b ,uj)F i (cr a ,(T 1 ,u)F^(a b ,a2,uj) 

^0 [fa, <7b, UJ J^oo Jqq Jqq 



(36) 



where 



i2a,ifc(x; <Ti,o- 2 ,w)= / d xi / d x 2 <( G o ,y(<7i,u;;xi)Go iW (0-2,u;;x 2 ; 
J© J© 



(5M ji (x 1 ,x 2 ,w) 



qo 



- / d 3 x" 
'0 



GoAmicTi, 00] X' 



Go(x",u;xi) Gnfcz(o-2,w;x 2 ) 



qo 



IIIJ 



+G ,ii(<Ti, w; x 1 )G* fcn (cr 2 , w; x A 



Go(x", u;; x 2 

qo / . 



M j/(xi,x 2 ,w) L (37) 



We note that the presence of three coupled volumetric integrals in eq. fl37j) would 
seem to make the calculation of the raw kernel computationally impractical. However, 
the functional derivatives of the linear wave operator generally will be proportional to 
three-dimensional spatial delta functions and their derivatives. Additionally, the reference 
source covariance matrix can be taken to be proportional to delta functions. Thus, the 
actual computational burden is much smaller than one might fear from the above result. 
We will see this in specific examples for kernels, such as those for sound speed and density. 



3. SOME SPECIFIC KERNEL FORMULAS 

In this section, we specify the aspects of acoustic wave propagation that we will 
consider and apply our general kernel formulas to a particular reference model of the Sun's 
structure and ensemble average of the source covariance matrix. We then give formulas for 
sensitivity kernels for perturbations in sound speed squared and density. 
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3.1. The Wave Equation and Reference Model 

In order to use the kernel formulas we have developed, we must make assumptions 
about the physics of acoustic wave propagation both in the Sun and in our reference model. 
Only then can we calculate our reference Green's functions and the functional derivatives 
of the wave equation operator that are needed. 

From here on, we will assume that there are no magnetic fields or bulk fluid flows 
present in the Sun. Additionally, we will work in the Cowling approximation, i.e., we 
assume that the density variations induced by the propagation of an acoustic wave have a 
negligible effect on the gravitational potential. For all but the lowest degree modes, this 
approximation is extremely accurate. In future work, we will relax most if not all of these 
restrictions. 



We will also i n clude damping in our wave equation. For simplicity, we will somewhat 



follow 



Birch et al. 



(120041 1 and use a simple convolutional model in the time domain. 
However, in contrast to their work, we will assume that the attenuation operator is spatially 
local. 



As for boundary conditions, we again follow iBirch et al.l ( 120041 ) and use a free surface 
boundar y condition, i .e., w e set the Lagrangian pressure equal to zero at some radius. As 



noted in 



Birch et al. 



( 120041 ). this does, of course, mean that waves with frequencies above 
the acoustic cutoff do not escape the Sun but are reflected back into it, in disagreement 
with reality. 



Under the above assumptions, the linearized equations for particle displacement u and 
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pressure perturbation P' have the following well-known form in the frequency domain: 

> , P 



-uTpu + VP' - g 



c 2 c 2 



u - (Vp) ■ u 



— iupTu = 0, 



P' + c 2 pV ■ u + pg • u = 0, 



(38) 
(39) 



where c and p are the sound speed and density, respectively, g is the gravitational 
acceleration, u is the angular frequency, and V is the damping parameter. For the real 
Sun, we assume that all of the medium parameters above can vary spatially in all three 
dimensions but not in time. In the case of the damping parameter, we assume that it also 
varies in frequency so that in the time domain it is convolved with the particle velocity v. 

However, for the purpose of constructing our wave equation operator £, we will rewrite 
eq. (J3HD slightly: 

1 _ cr 

u = 0. (40) 



- \p - V - u 



pCT Cr 

With this modification, our wave equation, eq. (125]) . becomes 



£ 







u 




P' 





iuT)pl - 4gg T + g(Vp) T V - 



_V T - — 



1 









u 




P' 



(41) 



where all vector quantities (and V) are taken to be column vectors and T indicates 
transposition. The reason for using eq. fj4*0|) is that in this form, the Green's matrix 
corresponding to £ satisfies reciprocity, i.e., 

G(x, w; x ) = G T (x , w; x), (42) 

Having a Green's matrix that obeys reciprocity greatly reduces the computational cost of 
computing kernels by limiting the number of source points for which we need to calculate 
responses, at least when one uses a technique other than mode summation, as we do. 

With a specific form of £ now given, we are in position to calculate its functional 
derivatives. In this paper, we will concern ourselves with kernels for sound speed squared 
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and density, so we only need the functional derivatives with respect to those quantities: 

<5£(x 



8c 2 (x') 



rS 



qo 



Pogog? 



T 

5o 



1 

Po. 



(43) 



AC fx) 



5p(x') 



5 fx 



^r )i - — gog? 



o 
i 



+ 



,V T <5 fx - x') 



(44) 



In eq. ( I44p we have neglected the derivative of g with respect to density under the 
assumption that the deviations in density from our reference model will cause such small 
changes in gravity that they will have a negligible effect on wave propagation. Note that 
this assumption is separate from the Cowling approximation. 

In addition to laying out our assumptions about wave propagation in the Sun, we must 
also describe our reference model. The obvious choice for a model of the Sun's structure 
is one that is spherically symmetric and in hydrostatic equilibrium. The symmetry of this 
choice greatly simplifies the calculation of the Green's functions we will need. 

As for the reference model of acoustic excitation, we will assume a simple form for 
the ensemb l e aver age of the source covariance matrix, eq. (13~T1) . We choose to generalize 



Birch et al. 



(120041 ) model to spherical coordinates: 

8(9 -9') 8 ((f) -4 



^My(x,x',u;)^ = S ir 8j 



-8' (r - r s ) 8' (r - r s ) f(u) 



(45) 



r 2 r' 2 sin 9 

where a subscript r indicates a component in the radial direction, and |/(a;)| 2 is the power 
spectrum of the excitation and may be adjusted to ensure that the waves modelled in the 
reference medium have an acoustic power spectrum close to the Sun's. Note that this a very 
practical model of excitation, not necessarily a realistic one. In particular, the presence of 
the ^-functions in the angular variables, which implies a zero horizontal correlation length 
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for the s ources, simplifie s the c omputational effort required to model the excitation. As 



noted in 



Gizon fc Birch 



(120021 ). the true correlation length of the sources is much smaller 



than the wavelengths of the waves we will model, anyway. 



3.2. Kernel Formulas for Perturbations to c 2 and p 

With our above assumptions in hand, we can now construct our sound speed and 
density kernel formulas. To begin, we introduce some simplifying notation. Recalling that 
our Green's function is a 4 x 4 matrix corresponding to a coupled system of equations in 
both particle displacement and pressure perturbation, we indicate its components by the 
indices r, 9, <p, and p, denoting the three directions in spherical coordinates and pressure. 
For instance, Go,rp(x, ou; xo) denotes the radial component of displacement at point x due 
to a point pressure source at x , in the reference medium, at frequency u. When we wish 
to leave a component unspecified (or indicate Einstein summation), we use Latin indices 
(other than r or p). 

To construct our kernel formulas, we start by incorporating the source excitation 
model, eq. (14"5]1 . into eq. (1371) . which we write as 

-R a ,;fc(x; °"1> °2, = Ha,ifc(x; 01, 02, ^) + #* jfei (x; CT 2 , (T 1 , id) (46) 

where 

Ha t ik(^] CT i,(T 2, Uj) = 



- / d 2 a'j <^jd'Go iim (tr u U)) x 

JdQ JO 



J 5C(x",u) 



5q a (x) 



d rs G (x.",uj;r s ,(T') ] d rs Gl M ((T 2 ,u;r s ,cr l ) 
qo 



(47) 

Note that d Ts is shorthand for partial differentiation with respect to r s , and we have 
assumed that the source covariance is fixed. 
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Now using eqs. f]4"3]) and (jHJ) and recognizing that go = —gor, we get 
F c 2(x; a 1 ,(T 2 ,uj) = 

--J dV 

c 0</<9© 

+ (5 ,r- P (o-i, w; x) ( -g- 5 r . s G'o ! rr(x, w; r s , cr') + — 9 rs Go,pr(x, u; r„ cr'] 

Po 



G 0|TT (o-i, w; x) ( p g%d rs Go !rr (x, w; r s , cr') - g d rs G 0:Pr (x, w; r„ cr r 



and 



# p (x; cr!, cr 2 ,u;) 



^2 / 

a cr 



90 



x9 r,Gn rr (ff2,w;r s ,«7') /(w 



+ iwr) ^ G ,rm(cri,a;;x)9 r . s G'o,mr(x, 



(4* 



a;;r s , cr 



Go, rr (c x)^|(9 rs (5o,rr(x, w; r s , cr') + G 0rp (cr X) oj] x) —^-^d Ts Go, pr (x, 



P6 c o 



u;r s ,cr 



+ <9 r ^Go,rr(cri, x)5f <9 r . s G'o,rr-(x, cu; r s , cr r 
+ -<9 e f G 0}rr (o"i, w; vL)g Q d r3 Go gr(x, u; r s , a') 



1 



( Go,rr(o"i, w ; x)5f 9 rs G ^r(x, w; r s , cr') 



r sin 6* 

x d rs G* 0rr (er 2 ,uj;r s ,cr') f(u) 
where all medium parameters are implicitly evaluated at x. 



(49) 



Lastly, we need to incorporate our source covariance model into our formula for D , 
the reference model version of eq. (129]) . which is needed in the evaluation of and C ' 
(eq. |22]) , which in turn are used in eq. ( 136]) . Using eq. ( 145]) , our formula for D becomes 

D (tr a ,a b ,u) = 



d a 1 d cr 2 / d (T , ,F i ((T a ,cr 1 ,u)d rs G 0jir {(T 1 ,uj]r s ,(T') 

90 Jd& JdG 

2 

x F^{a h ,a 2 ,u))d rs Gt kr (<r 2 ,uj; r s ,a') f(u) . (50) 



Thus, to construct traveltime sensitivity kernels for squared sound speed and density, 
under the assumptions laid out above, the equations that are needed are ( 120]) . (122"]) . (136]) . 
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( 146]) . and f l48|) - fl50|) . Of course, the user of these results will still need to provide additional 
information, such as a specific reference model, observational geometry, and so forth. 



EXAMPLES 



In this section we display examples of kernels for squared sound speed for large 
distances, in order to give an idea of the extent and character of traveltime sensitivity 
at depth in the Sun. In order to accelerate the calculations somewhat, we make one 
particularly simplifying assumption, i.e., we take our filter function F (see eq. [6]) to be 



F(tr,(T , ,t) = f5(tr-tT , )5(t) 



(51) 



The assumption that the line of sight is always radial is quite helpful in lessening the 
computational complexity, but it is, of course, only a good approximation near the center 
of the Sun's visible disk. 



As for our reference model, we use model S (jChristensen-Dalsgaard et al.lll996l ) as the 
model of the Sun's structure, and a free-surface boundary condition, as mentioned earlier. 
For the power spectrum of the acoustic excitation, we once again turn to 
and use 



Birch et al. 



feooj) 



/(") 



-w 2 T 2 



(52) 



with To, 



68 s. 



Our damping model is actually a fusion of the ones used in 



Gizon &: Birch 



Birch et al. 



(120041 ) . We let 



T(u,r) = g(u)h(r), 



where from 



Gizon fc Birch 



(120021 ) we take 



g(u) = 7 



(J200J) and 

(53) 

(54) 
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with 7/27T = 100/iHz, u*/2it = 3mHz, and /3 = 4.4. Whereas from iBirch et all (120041 ) we 

get 

r < T M-^ 2 1 (55) 



/i(r) = exp 
with 



(At) 



2 



r dr' 

T(r) = / — -, (56) 

Jo c o(n 

T c = T(R), where R is the radius of the photosphere, and At =69 s. 

Lastly, we assume that our time window W is simply one for times within 10 m of a 
first arrival and zero otherwise. 

In Fig. 1, we see our first example. Shown is a kernel for squared sound speed for two 
observation points separated by and angle of 30°. Because of the huge range of amplitudes 
in the kernel, the strongest regions being close to the observation points, we have scaled the 
kernel by the sound speed and saturated the scale in order to bring out the deep structure. 
Indicated in both cross-sections is the location of the ray path associated with the first 
arrival at this distance. Note that the kernel displays the classic "banana-doughnut" shape, 
with sensitivity reaching a local minimum on the ray path. The images here are limited to 
the upper 40% of the Sun, by radius. We can see that the kernel has sensitivity mostly in 
the convection zone and above. 

In Fig. 2, we have the squared sound speed kernel for a distance of 60°. In this case, 
the sections show the upper 70% of the Sun. Clearly, the kernel in this case significantly 
intersects the tachocline, indicating that it may be of some use in interpreting time-distance 
measurements aimed at studying that region. 

In both figures, it is clear that ray theory cannot be a good approximation here, given 
that the waves are clearly sensitive to a broad volume around the ray. 
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Fig. 1. — Vertical cross-sections of a sensitivity kernel for squared sound speed for a distance 
of 30°. The kernel has been scaled by the sound speed and the scale has been severely 
saturated. (Top) The cross-section in the plane of the ray path (black line). (Bottom) The 
cross-section perpendicular to the ray path. The intersection point of the plane and the ray 
is indicated by the black dot. The total angular range in both sections is 50°, and both 
extend radially from 0.6 R & to the surface. Note that this is the kernel for waves travelling 
from point B to point A. The units are Mm~ 2 
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Fig. 2. — Vertical cross-sections of a sensitivity kernel for squared sound speed for a distance 
of 60°. Again the kernel has been scaled by the sound speed and the scale has been saturated. 
( Top) The cross-section in the plane of the ray path (black line). (Bottom) The cross-section 
perpendicular to the ray path. The intersection point of the plane and the ray is indicated 
by the black dot. The total angular range in both sections is 100°, and both extend radially 
from 0.3 R Q to the surface. 



5. DISCUSSION 



We have presented a different theoretical framework, based on functional analytic tools, 
that we hope many will find a clear and elegant way of thinking about traveltime sensitivity 
kernels. In addition, we have provided the essential formulas for calculating these kernels in 
a full spherical geometry. For kernels for medium properties other than the ones discussed 
here, the mathematical apparatus has been developed to allow the reader to derive his or 
her own formulas. 

With some reasonable assumptions about the Sun, we have shown that one can obtain 
kernels that demonstrate their potential value by showing the inadequacy of ray theoretical 
approximations. Clearly these results can now be put to good use in analyzing time-distance 
measurements to obtain information on the deep structure of the Sun. 
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